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SUMMARY 


A correlation procedure which relates strength to flaw 
'’magnitude'’ is reviewed and presented for ceramic mater- 
ials. The procedure was applied to 20 nominally identical 
polished polycrystalline aluminum- oxide flexure specimens. 
Using the maximum linear dimension of pores and grain pull- 
outs in the flexure specimens as a measure of the stress- 
concentrating ability of surface flaws, an exponential distri- 
bution was fitted to an observed cumulative flaw "magnitude” 
distribution, and its associated largest value distribution 
calculated. These calculated values were compared with 
values which were observed on gauge-size pieces of failed 
samples. Using this largest value distribution, flaw mag- 
nitude and strength were correlated and the results compared 
with the Griffith equation. The comparison showed that equa- 
tions with the Griffith form describe the data reasonably 
well. The comparison served to define clearly the correla- 
tion procedure and problems that are associated with its 
application. 
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PHENOMENOLOGICAL AND STATISTICAL ANALYSIS OF 
FRACTURE IN POLYCRYSTALLINE ALUMINUM OXIDE 


by 

Thomas M. Heslin and Alfred G. Eubanks 
Goddard Space Flight Center 


INTRODUCTION 

According to the ’’weakest -link” model, a material has a unique strength called its ’’theoretical 
strength,” which is pfthe order of 10® psi for ceramics. This strength does not in general actually 
exist, because of the presence of numerous flaws, each of which acts as a stress concentrator. 
Fracture occurs when the peak stress adjacent to the single, most acute flaw in a body reaches the 
’’theoretical strength” of the material. Of course, this strength is purely conceptual in nature, as 
one cannot test the strength of an isolated element without changing the conditions that exist when 
the element is actually in the body. In addition to this indeterminacy, there is an uncertainty as- 
sociated with the statistical nature of the concept itself which is considered by some investigators 
to be an inherent property of materials. This inherent uncertainty arises because the acuteness of 
the most severe flaw in a specimen may differ for individual specimens in a nominally identical set. 
Thus, if flaw ’’magnitudes” are statistically distributed in the parent body from which a set of test 
samples is taken, there is a second statistical distribution, related to the first, for the most severe 
flaw which exists in samples taken from the parent (Reference 1). It is generally assumed that this 
distribution is directly related to an observable relationship commonly referred to as the fracture 
probability. The fracture probability, denoted S(a), is defined (Reference 2) by 

’ ( 1 ) 


where N is the total number of samples tested; m is the specimen serial number corresponding to 
a list of fracture stresses arranged in an increasing order from 1 to N; a is the observed fracture 
stress for a given m. 

The mathematics involved in the weakest- link model is identical to the sampling problem in 
statistics of finding the least value in samples of size n, drawn from a population with a known 
value distribution (Reference 3). 

If the stress-concentrating ability, or ’’magnitude,” of a flaw is denoted by ”c”, and strength 
is assumed to be an inverse function of c, then strength and flaw ’’magnitude” can be related by 


finding a cumulative distribution, F(l/c), for the flaw population in the parent body. Given F(l/c), 
a cumulative distribution of the smallest value, G(l/c) , in samples of size n is given (Reference 3) 
by 


G 




( 2 ) 


It is assumed that this distribution is directly related to the observed fracture probability, 
S(o-), and that strength and flaw ^’magnitude^’ can be correlated by comparing i/c and a- for equal 
values of G(l/c) and S(cr). 

The study objectives were to establish the types of surface flaws present in a polycrystalline 
aluminum oxide compact by thoroughly characterizing samples and to define more clearly the pro- 
cedure outlined and the problems associated with its application. 


CHARACTERIZATION AND TESTING 

Twenty 5” by 1/2” by 1/8” flexure specimens of 99.5 percent AI 2 O 3 were procured. A tabu- 
lation of various characterization parameters is given in Table 1. The 5” by 1/2” surfaces of the 


Table 1 


Quantitative Characterization of Flexure Specimens. 


Specimen 

Serial 

Number 

Modulus of 
Rupture x 10“^ 
(psi) 

Porosity 

{%) 

Average 
Grain Size 
W 

Average Linear 
Flaw Density 
(counts/ 
inch) 

Average 

(CLA) 

(micro- 

inches) 

Average 
Areal Flaw 
Density x 10"^ 
(counts/ 
sq. in.) 

Largest 
Flaw X 10^ 
(inches) 

1 

26.959 

4.378 

24.43 

256 

4.27 

9.54 

4.43 

2 

27.807 

4.128 

23.67 

392 

6.69 

12.9 

5.00 

3 

27.952 

4.413 

23.67 

357 

7.11 

10.3 

6.25 

4 

29.022 

3.840 

16.11 

251 

3.17 

10.4 

5.92 

5 

29.563 

5.477 

17.85 

388 

4.87 

11.0 

5.06 

6 

29.720 

4.242 

25.53 

237 

3.41 

10.9 

5.16 

7 

30.012 

4.245 

20.32 

194 

3.91 

8.11 

5.58 

8 

30.452 

4.213 

26.56 

252 

4.21 

9.05 

5.21 

9 

31.078 

4.249 

24.61 ' 

299 

5.31 

6.46 

6.00 

10 

31.168 

4.149 

17.80 

281 

5.20 

5.78 

4.95 

11 

31.315 

4.357 

13.13 

287 

5.00 

7.12 

1 5.93 

12 

31.908 

4.038 ! 

23.02 

408 

5.18 

18.8 

5.81 

13 

32.937 

4.459 1 

26.87 

240 

1 4.26 

1 9.09 

5.16 

14 

34.675 

5.037 ' 

11.05 

439 

, 5.19 

11.1 

5.09 

15 

34.721 

, 5.223 

11.03 

397 

3.78 

11.1 

5.46 

16 

34.846 

5.536 

11.73 

' 448 

4.29 

12.6 

4.62 

17 

34.914 

i 5.011 

11.99 

362 

4.72 

10.7 

5.10 

18 

i 34.948 

5.528 

14.96 

455 

4.20 

13.8 

5.17 

19 

36.181 

4.348 

j 15.60 

221 

3.06 

6.89 

5.77 

20 

36.314 

5.571 

12.34 

381 

4.31 

12.8 

4.86 
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samples were polished by the vendor and had surface roughnesses ranging from 3.06 microinches 
to 7.11 microinches CLA (center line average). Photomicrographs of a sample, as received and 
after thermal etching, are shown in Figure 1. The samples were failed in the four-point flexure 
fixture shown in Figure 2. This fixture is designed to meet the specifications of ASTM method 
C78-64. The gauge region is 1-1/3” long. The fixture is activated by placing it in the compression 
stroke of an Instron Model TTC testing machine. The top and bottom of the fixture are attached 
magnetically to the Instron machine's traveling head and compression cell, respectively. The side 



Before thermal etching 
21 5x 



21 5x 

Figure 1— -Sample 15 before and after thermal etching. 



Figure 2—FaiIure fixture. 


posts serve only to initially align the fixture, 
and are dropped free of the action before 
testing. 

A spectrochemical analysis for mag- 
nesium and chromium additives of the as- 
received material was performed and showed 
chromium and magnesium present in aver- 
age amounts of 0.23 wt percent and 0.17 wt 
percent, respectively. Bulk porosities of 
the samples were determined by comparison 
of measured density to the theoretical den- 
sity of 3.987g./cc. The porosity ranged from 
a high of 5.571 percent to a low of 3.840 per- 
cent. Grain size was measured by thermally 
etching pieces of failed specimens, photo- 
graphing them, and measuring the intercepts 
of 120 grains with lines drawn across the 
photomicrograph. Grain size ranged from a 
maximum of 26.87 microns to a minimum of 
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11,03 microns. Areal flaw density was measured by counting the number of pores and grain pull- 
outs in two 215 X photomicrographs of each specimen. One of the photographs was taJcen near the 
center of the gauge area; the other was taken near the long edge of the gauge area on the tension 
stirface. The flaw densities ranged from a maximum of 1.88 x 10® flaws/in.^ to a minimum of 
5.78 X 10® flaws/in. ^ Linear flaw density was measured by means of a Talysurf 4 profilometer 
and by counting the spike-like depressions in two profiles, each of which had 10,000 x vertical mag- 
nification, 500 X horizontal magnification, and a 0.3-inch stylus traverse. A representative profile 
example is shown in Figure 3. One profile was taken near the center of the gauge area, while the 
other was taken near the long edge of the area. Linear flaw density ranged from a maximum of 455 
flaws/in. to a minimum of 221 flaws/in. 

A piece of a failed specimen was submitted to X-ray analysis to determine whether there was 
any preferred orientation of surface crystallites. None was found. Fracture surfaces and areas 
adjacent to fracture surfaces were examined with a scanning electron microscope (SEM); repre- 
sentative pictures are shown in Figures 4 and 5. Replicas were made of fracture surfaces and ex- 
amined with a transmission electron microscope. Representative photographs of the surface using 
this technique are shown in Figure 6. Some samples had hairline cracks branching off from the 
main fracture. One such crack was examined before and after thermal etching, and is shown in 
Figure 7. Another such crack was examined with the scanning electron microscope before and 
after thermal etching; it is shown in Figure 8. Pits at points of dislocation termination were ob- 
served on individual grains after thermal etching. An example of this type of flaw is shown in 

10,000x vertical magnification 
500x horizontal magnification 
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Before thermai etching 
195x 



After thermal etching 
195x 


Figure 7— Hairline crack before and 
after thermal etching. 



Before thermal etching 
900x 



After thermai etching 
900x 


Figure 8— Hairline crack before and after thermal etching 
viewed through SEM. 


Figure 9. The square geometry of some of the pits indicates that magnesium is present as mag- 
nesia and that there are grains of the magnesia dispersed in the body. 

The above characterization techniques revealed the presence of five types of flaws in these 
bodies. They are: pores; grain pullouts from polishing; polishing scratches; grain boundaries; 
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lOOOx 


'■»» V*' si f 


“K , 


^ 

due to pores and grain pullouts. The frac- 
lOOOx & ^ 

ture path in failed specimens seemed to 

~ ^ include as many of these holes as possible. 

X / I Grain boundaries also are a pronounced 

i ^ feature, but only after the polished damage 

? “is \-' " f^****! layer has been removed by thermal etching. 

4 "" / The damage layer probably has tensile 

properties that are quite different from the 
material; and it is difficult to deter- 
\ mine whether fracture actually nucleates in 

j this layer, or just below it at the interface 

? ^ grain boundary and the damage layer. 

Grain boundaries are a prime candidate 

I , consideration as fracture-nucleating 

'X' flaws. It has been fairly well established 

that thermal expansion and elastic aniso- 
lOOOx ^ 

tropy at grain boundaries are major causes 

Figure 9-Pictures showing etch pit in individual grains. weakness in macroscopic bicrystals (Ref- 

erences 5 and 6). However, on the micro- 
scopic scale of a poly crystalline compact, the forces across the curved peripheral boundaries be- 
tween particles, which drive the sintering process and tend to pull the grains together, would be of 
great significance; but such forces would be relatively unimportant for a macroscopic bicrystal. 

The cohesive force system between particles in the two instances is somewhat different. There is 
thus some uncertainty in directly extrapolating results obtained with macroscopic samples to a 
microscopic situation. 

It was impossible to firmly classify fracture of these 20 samples as occurring between or 
across grains, because instances of both types of propagation were observed. Considering these 
facts, and for lack of a more prominent irregularity on unetched, polished surfaces, it was assumed 


Figure 9— Pictures showing etch pit In Individual grains. 
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that pores and grain pullouts were fracture -nucleating flaws. The most obvious and easily made 
measurement of pores and grain pullouts is their maximum linear dimension as measured from 
photomicrographs. It was assumed that this dimension is an indication of the stress -concentrating 
ability of pores and grain pullouts. 

The flaw measurement procedure was as follows. Two 215 x photomicrographs were taken in 
the gauge area of each of the 20 samples. A grid was then put over the pictures, and the maximmn 
linear dimension of each pore and grain pullout was measured using a pair of dial calipers. These 
measurements ranged from 0.01 to 0.85 inches. The measurements were arranged in increasing 
order; and the number of flaws, n, in each 0.010 inch interval was taken as the number of flaws 
with magnitude equal to the center value of that particular interval. Flaw magnitude values, de- 
noted c, were calculated, and values of l/c tabulated along with the cumulative number, 5:n, of flaws 
with value less than l/c. This information is shown in Table 2. A cumulative flaw distribution was 
found by dividing values of Xn by the total number of flaws, N, counted in all the photographs. This 
cumulative flaw distribution was taken as representing the parent flaw distribution, F( l/c) . 


Table 2 


Tabulation for a Cumulative Flaw Distribution for 20 Specimens. 


rval Center 

c X 10"* 

7 X 10"2 


>'n 

>'n 

{ >n\ ^ 

(in.) 

(in.) 

(in.'i) 

n 

N 

‘°8 --N/ 

0.015 

0.698 

143 

1428 

14,195 

1.00 

'.X) 

0.025 

1.16 

86.0 

2429 

12,767 

0.899 

0.997 

0.035 

1.63 

d1.4 

2109 

10,338 

0.728 

0.566 

0.045 

2.09 

47.8 

1775 

8229 

0.580 

0.376 

0.055 

2.56 

39.1 

1408 

6454 

0.455 

0.263 

0.065 

3.02 

33.1 

1051 

5046 

0.355 

0.190 

0.075 

3.48 

28.7 

820 

3995 

0.281 

0.144 

0.085 

3.95 

25.3 

592 

3175 

0.224 

0.109 

0.095 

4.42 

22.6 

443 

2583 

0.182 

8.71 X 10'^ 

0.105 

4.88 

20.5 

336 

2140 

0.151 

6.08 X 10 ‘2 

0.115 

5.35 

18.7 

268 

1804 

0.127 

5.88 X 10'^ 

0.125 

5.81 

17.2 

208 

1536 

0.108 

5.00 X 10'^ 

0.135 

6.28 

15.9 

154 

1328 

9.35 X 10"^ 

4.26 X 10 

0.145 

6.74 

14.8 

151 

1174 

8.27 X 10'^ 

3.74 X 10’^ 

0.155 

7.21 

13.9 

112 

1023 

7.21 X 10‘^ 

3.26 X 10“^ 

0.165 

7.67 

13.0 

100 

911 

6.42 X 10’^ 

2.86 X 10"^ 

0.175 

8.14 

12.3 

90 

811 

5.71 X 10‘^ 

2.53 X 10“^ 

0.185 

8.60 

11.6 

71 

721 

5.08 X 10'^ 

2.24 X 10' 2 

0.195 

9.07 

11.0 

63 

650 

4.58 X 10'^ 

2.04 X 10"^ 

0.205 

9.53 

10.5 

60 

587 

4.13 X 10 

1.83 X 10'2 

0.215 

10.0 

10.0 

47 

527 

3.71 X 10‘2 

1.62 X 10’^ 

0.225 

10.5 

9.55 

39 

480 

3.38 X 10'^ 

1.49 X 10'^ 

0.235 

10.9 

9.15 

28 

441 

3.11 X 10"^ 

1.37 X 10' 2 

0.245 

11.4 

8.77 

34 

413 

2.91 X 10’^ 

1.26 X 10'^ 

0.255 

11.9 

8.43 

25 

379 

2.67 X 10'^ 

1.16 X 10'^ 

0.265 

12.3 

8.11 

27 

354 

2.49 X 10'^ 

1.08 X 10'^ 

0.275 

12.8 

7.82 

22 

327 

2.30 X 10‘2 

1.00 X 10'2 

0.285 

13.3 

7.54 

23 

305 

2.15 X 10’^ 

9.33 X 10'^ 
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Table 2 (Continued) 


Interval Center 

c X 10“* 

i X 10-2 


2n 

2n 

/ 2n\‘‘ 

(in.) 

(in.) 

(in.-‘) 

n 

N 

log (i - -pr) 

0.295 

13.7 

7.29 

21 

282 

1.99 X lO'^ 

8.64 X 10-3 

0.305 

14.2 

7.05 

21 

261 

1.84 X 10-2 

7.98 X 10-3 

0.315 

14.6 

6.82 

16 

240 

1.69 X 10-2 

7.33 X 10-3 

0.325 

15.1 

6.61 

21 

174 

1.22 X 10-2 

5.29 X 10-3 

0.335 

15.6 

6.42 

12 

153 

1.07 X 10-2 

4.64 X 10-3 

0.345 

16.0 

6.23 

23 

133 

9.37 X 10-^ 

4.07 X 10-3 

0.355 

16.5 

6.05 

15 

110 

7.75 X 10 

3.36 X 10-3 

0.365 

17.0 

5.89 

8 

95 

6.69 X 10-^ 

2.90 X 10-3 

0.375 

17.4 

5.73 

9 

86 

6.06 X 10"^ 

2.63 X 10-3 

0.385 

17.9 

5.58 

10 

77 

5.42 X 10" 2 

2.35 X 10"3 

0.395 

18.4 

5.44 

6 

67 

4.72 X 10"^ 

2,04 X 10"^ 

0.405 

18.8 

5.31 

7 

61 

4.30 X 10'^ 

1.87 X 10"3 

0.415 

19.3 

5.18 

8 

54 

3.80 X 10"^ 

1.65 X 10-3 

0.425 

19.8 

5.06 

6 

46 

3.24 X 10"2 

1.41 X 10-3 

0.435 

20.2 

4.94 

4 

40 

2.82 X 10"^ 

1.23 X 10"3 

0.445 

20.7 

4.83 

4 

36 

2.53 X 10"^ 

1.10 X 10-3 

0.455 

21.2 

4.72 

3 

32 

2.25 X 10"^ 

9.75 X 10"“ 

0.465 

21.6 

4.62 

5 

29 

2.04 X 10-=* 

8.85 X 10"“ 

0.475 

22.1 

4.53 





0.485 

22.5 

4.43 

3 

24 

1.69 X 10"2 

7.33 X 10"“ 

0.495 

23.0 

4.34 





0.505 

23.5 

4.26 

1 

21 

1.48 X 10"^ 

6.41 X 10-“ 

0.515 

24.0 

4.17 

2 

20 

1.41 X 10"^ 

6.11 X 10"“ 

0.525 

24.4 

4.09 

2 

18 

1.27 X 10"3 

5.51 X 10-“ 

0.535 

24.9 

4.02 





0.545 

25.3 

3.94 





0.555 

25.8 

3.87 





0.565 

26.3 

3.80 

1 

16 

1.13 X 10"2 

4.90 X 10-“ 

0.575 

26.7 

3.74 





0.585 

27.2 

3.67 

2 

15 

1.06 X 10"2 

4.60 X 10"“ 

0.595 

27.6 

3.61 

2 

13 

9.16 X lO"'* 

3.98 X 10-“ 

0.605 

28.1 

3.55 





0.615 

28.6 

3.50 

1 

' 11 

7.75 X 10 

3.36 X 10"“ 

0.625 

29.1 

3.44 

1 

10 

7.04 X 10"'' 

3.06 X 10-“ 

0.635 

29.5 

3.39 





0.645 

30.0 

3.33 

2 

9 

6.34 X lO-"* 

2.75 X 10-“ 

0.655 

30.5 

3.28 

1 

7 

4.93 X 10" 

2.14 X 10'“ 

0.665 

30.9 

3.23 

1 

6 

4.23 X 10" 

1.84 X 10-“ 

0.675 

31.4 

3.18 





0.685 

31.9 

3.14 





0.695 

32.3 

3.09 





0.705 

32.8 

3.05 





0.715 

33.2 

3.01 





0.725 

33.7 

2.96 

1 

5 

3.52 X lO"* 

1.53 X 10"“ 

0.735 

34.2 

2.92 

1 

4 

2.82 X 10-'* 

1.22 X 10"“ 

0.745 

34.6 

2.89 





0.755 

35.1 

2.85 

1 

3 

2.11 X 10-* 

9.15 X 10"* 

0.765 

35.6 

2.81 

1 

2 

1.41 X 10-“ 

6.11 X 10-3 

0.775 

36.0 

2.77 





0.785 

36.5 

2.74 





0.795 

37.0 

2.70 





0.845 

39.3 

2.54 

1 

1 

7.04 X 10-® 

3.05 X 10" 3 
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The average flaw density in the samples was thus calculated to be 1.04 x 10® flaws per inch^. 
This gives a value of 6.92 x 10® flaws per gauge area. Thus, for G(l/c) = 0.01, F(l/c) is found by 
Equation 2 to be 1.59 x 10“®, and for G(l/c) = 0.99, F(l/c) is 6.65 x 10”®. Therefore, in order to 
determine a G(l/c) distribution without having to extrapolate F(l/c) , experimental values of F(l/c) 
must be found on the order of 10”^. This means that a total flaw count on the order of 10® or 10® 
flaws should be made. To develop a rapid method of measuring the stress-concentrating ability of 
a flaw with flaw densities in the ranges observed is a problem, since, at a measuring rate of one 
flaw per second, it would require on the order of 3 to 30 years to record the necessary data. 

As the desired experimental values of F(l/c) could not be found, the data which were obtainable 
in a reasonable amount of time were fitted to a number of exponential distributions. The simplest 
one, and also one that gave a good fit, was 

f(-^) = 1 - = , (3) 

where a and are parameters determined by plotting log log [l - (2n/N)]"^ vs log (l/c). A tabula- 
tion of l/c and log [l - (Sn/N)] is given inrTable 2, and a plot of log log [l - (Xn/N)]"^ vs log ( l/c) 
is shown in Figure 10. For small values of Sn/N, the distribution 1 - with a = 2.20 x 10“^® 

and 13 = 4.88 was found to represent F(l/c) well, and was used to extrapolate values of l/c for the 
distribution of smallest values, G( l/c) , using Equation 3 in Equation 2. The results are shown in 
Table 3. ! ’• 

Table 3 


Extrapolation of c Values for G(l/c) Using Equation 3 in Equation 2. 




1 

/ 1 

1 (\\^ 

1 


Glci 

- In(l-G) 

X In(l-G) 


2. 3/3 X 

c 

C X 10® 

0,048 

0.04919 

3.23 X 10® 

19.59312 

1.74144 

55.14 

18.1 

0.095 

0.09982 

6.56 X 10® 

20.30163 

1.80441 

63.69 

15.7 

0.143 

0.15432 

1.01 X 10® 

20.73317 

1.84276 

69.62 

14.4 

0.190 

0.21072 

1.38 X 10® 

21.04530 

1.87051 

74.22 

13.5 

0.238 

0.27181 

1.78 X 10® 

21.29983 

1.89313 

78.19 

12.8 

0.286 

0.33687 

2.21 X 10® 

21.51621 

1.91236 

81.73 

12.2 

0.334 

0.40647 

2.67 X 10® 

21.70530 

1.92917 

84.95 

11.8 

0.381 

0.47965 

3.15 x,10® 

21.87062 

1.94386 

87.87 

11.4 

0.428 

0.55862 

3.67 X 10® 

22.02341 

1.95744 

90.56 

11.0 

0.476 

0.64626 

4.25 X 10® 

22.16778 

1.97027 

93.38 

10.7 

0.523 ! 

0.74024 1 

4.86 X 10® 

22.30426 

1.98240 

96.03 

10.4 

0.570 

0.84397 

5.55 X 10® 

22.43521 

1.99404 

98.64 , 

10.1 

0.618 

0.96233 

6.32X10® 

22.56694 

2.00575 

101.3 

9.87 

0.666 

1.09362 1 

7.18x10® 

22.69450 

2.01709 

104.0 : 

9.61 

0.714 

1.25176 

8.22 X 10® 

22.83101 

2.02922 

107.0 

9.34 

0.716 

1.43129 

9.40 X 10® 

22.96499 ; 

2.04113 

109.9 

9.10 

0.810 

1.66073 

1.09 X 10^® 

23.11198 

2.05488 

113.5 

8.81 

0.850 

1.89712 

1.25 X lO^® 

23.24092 

2.06565 

116.3 

8.60 

0.900 j 

2.30259 , 

1.51 X 10‘° 

23.43791 

2.08316 

122.1 

8.19 

0.952 1 

3.63655 

1.99 X 10 

23.71394 

2.10769 

128.1 

7.81 


Average number of flaws per gauge area = 6.92 x 10® 
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The samples were failed and, using Equation 1, the fracture probability, S(o-), was computed. 
A functional relation between a and l/c of the form cr = m( 1/c)p was assumed, and for equal values 
of S(cr) and G(l/c) the value of log (a) was plotted vs log (l/c). The result is shown in Figure 11. 
A tabulation of log (o-) and log (l/c) for equal values of S(o-) and G(l/c) is given in Table 4. 

Gauge- size pieces of failed specimens were scanned with a microscope, and the maximum 
linear dimension of the largest pore or grain pullout was recorded. These values are given in 
Table 1. 



Figure 11— Plot of log o* vs Iog(l/c) for equal values of S(ct) and G(l/c). 
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Table 4 


Tabulation of log(c^ and log(l/c) for Equal Values of S(a*) and G(l/c). 


S(a) or G(l/c) 

o-x 10"^ 

log(cr) 

1/c 

log(l/c) 

0.048 

26.959 

4.43249 

55.14 

1.74144 

0.095 

27.807 

4.44514 

63.69 

1.80441 

0.143 

27.952 

4.45576 

69.62 

1.84276 

0.190 

29.022 

4.46389 

74.22 

1.87051 

0.238 

29.563 

4.47070 

78.19 

1.89313 

0.286 

29.720 

4.47712 

81.73 

1.91236 

0.334 

30.012 

4.48287 

84.95 

1.92917 

0.381 

30.452 

4.48855 

87.87 

1.94386 

0.428 

31.078 

4.49248 

90.56 

1.95744 

0.476 

31.168 

4.49776 

93.38 

1.97027 

0.523 

31.315 

4.50365 

96.03 

1.98240 

0.570 

31.908 

4.50840 

98.64 

1.99404 

0.618 

32.937 

4.51388 

101.3 

2.00575 

0.666 

34.675 

4.51943 

104.0 

2.01709 

0.714 

34.721 

4.52530 

107.0 

2.02922 

0.761 

34.846 

4.53212 

109.9 

2.04113 

0.810 

34.914 

4.53857 

113.5 

2.05488 

0.850 

34.948 

4.54507 

116.3 

2.06565 

0.900 

36.181 

4.55328 

122.1 

2.08316 

0.952 

36.314 

4.56086 

128.1 

2.10769 


Assume cr - M(1/c)p. 


Then log cr = p log ( 1/c) + log M. 


RESULTS AND DISCUSSION 


The Griffith equation for fracture strength (Reference 7) is 


a 


/ 


77 Ey 

cir^) 


where 

E = modulus of elasticity, 

y = specific surface energy of new surfaces formed, 
c = major diameter of an elliptical crack, and 
V = Poisson^ s ratio. 

Values of E and v are given (Reference 8) as 52 x 10® psi and 0.21, respectively. The surface 
energy, y, associated with the fracture process generally consists of a number of components among 
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which is the chemical surface free energy. Recent values of y reported for single crystalline 
material range from 6.0 to 7.3 J/m^ depending on crystal orientation (Reference 9). Using an 
average value of 6.6 J/m^ or 3.76 X in.-lbs/in.^, one gets cr = 2.54 x 10^ (l/c)^/^ ^ 

expressed in inches and o- expressed in pounds per square inch. 

A rough estimate of the applicability of an equation of the Griffith form using this analysis is 
offered in Figure 11. As shown, the data depict more of a curve than a straight line; however, the 
curve is well defined by two lines. The upper line has a slope and an intercept which are in order 
with those of the Griffith equation. 

If the strength of a material is governed by a surface flaw distribution, and if one considers 
only the uniformly stressed central portion of a beam under four -point loading, then the Weibull 
equation for probability of fracture (Reference 10) is 


where 


A = the gauge area of samples tested, 

= the actual fracture stress, 
o-^ = the stress below which fracture cannot occur, 

= a material constant, and 
m = a constant representative of the flaw density. 

The Weibull equation with = 0 and Equation 3 are quite similar if a is substituted for l/c 
using an equation of the Griffith form. The result is 

.t m ^ 

F(cr) = 1 - e'A , 


where A' = a/M^/p and m' ^ /3/p. 

The shape of the curve in Figure 10 is interesting because of the kink in it. This suggests a 
bimodal distribution of flaws. Pores and grain pullouts are visually indistinguishable in most 
cases, and no attempt at distinguishing between them was made while counting flaws for this anal- 
ysis. It is therefore possible that this kink is the result of superposing the two flaw distributions. 

CONCLUSIONS 

A thorough characterization of 20 high-alumina (99.5 percent AI 2 O 3 ) flexure specimens es- 
tablished the presence of the following types of flaws: 

grain pullouts due to polishing 
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polishing scratches 

pores 

grain boundaries 

dislocations in individual grains. 

Application of a statistical correlation procedure which relates strength to flaw "magnitude" 
demonstrated the need for considerable extrapolation. The fact that there is order -of -magnitude 
agreement between observed and calculated values of c(max) indicates an internal consistency, and 
therefore the assumptions made were reasonably accurate. Considering the clarity with which 
S(cr) was defined using only 20 specimens, the fact that a plot of log (or) vs log(i/c) for equal values 
of S(a) and G(l/c) is a continuous curve and compares with the Griffith equation indicates that a 
more exhaustive investigation which systematically attacks some problems brought out by this in- 
vestigation may prove frmtful (at least in the area of material predictability). 

This investigation has established the following problems: 

1. Determining what type of flaw (if indeed there is just one type) nucleates fracture in a 
brittle poly crystalline material, such as a sintered AI 2 O 3 compact, is a problem. Un- 
doubtedly there is interaction between flaw networks, and such interaction could preclude 
separation of the networks and prevent superposing their effects. 

2. Determining a measurement or series of measurements which significantly characterizes 
the stress- concentrating ability of existing flaws is a problem. 

3. Devising a counting device which measures the "magnitude" of a flaw and does it at such a 
rate as to make practical an experimental determination of parent flaw distributions is a 
problem. 

Solution of the last of these problems is necessary to the confident application of a statistical 
correlation procedure based on the "weakest-link" concept. Actually, if a parent flaw population 
could be adequately characterized experimentally, such a procedure could be used to pragmatically 
test possible solutions to the first two problems. This, of course, assumes applicability of a 
weakest-link type theory in which, conceptually at least, one is dealing with the breakage of a chain 
with links joined in series with the force on each link equal to the force applied to the chain as a 
whole. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, December 1 , 1969 
601 - 13 - 02 - 42-51 
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